Microscopic theory of the proximity effect in superconductor-graphene nanostructures 
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We present a theoretical analysis of the proximity effect at a graphene-superconductor interface. 
We use a tight-binding model for the electronic states in this system which allows to describe the 
interface at the microscopic level. Two different interface models are proposed: one in which the 
superconductor induces a finite pairing in the graphene regions underneath, thus maintaining the 
honeycomb structure at the interface and one that assumes that the graphene layer is directly coupled 
to a bulk superconducting electrode. We show that properties like the Andreev reflection probability 
and its channel decomposition depend critically on the model used to describe the interface. We 
also study the proximity effect on the local density of states on the graphene. For finite layers we 
analyze the induced minigap and how it is reduced when the length of the layer increases. Results 
for the local density of states profiles for finite and semi-infinite layers are presented. 

PACS numbers: 73.23.-b, 74.45. +c, 74.78.Na, 73.20.-r 
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INTRODUCTION 

The possibility to isolate and perform direct transport 
measurements on few or even single graphite layers [l| 
has triggered a large activity in the condensed matter 
community. The case of a single layer of carbon atoms, 
known as graphene, is of particular interest because of 
its unique electronic structure which, under certain con- 
ditions corresponds to massless Dirac fermions confined 
in two dimensions 0. 

On the other hand, the coupling to a superconductor 
provides an interesting way to test the electronic prop- 
erties of graphene. In a recent work by Beenakker 0] it 
was shown that for an ideal interface between a supercon- 
ductor and graphene an unusual type of Andreev reflec- 
tion, in which the hole is specularly reflected, appears. 
Several other effects involving graphene and supercon- 
ductors like Josephson transport 0, Q , re-entrance effect 
Q , and quasiparticle transport mediated by multiple An- 
dreev processes |7[] have been theoretically analyzed. 

In addition to its effect on the transport properties, 
the coupling to a superconductor also should produce a 
change in the electronic spectral properties and the in- 
duction of pairing correlations due to the proximity ef- 
fect. The recent experimental achievement of good con- 
tact between superconducting electrodes and graphene 
layers [|| open the possibility to explore the proxim- 
ity effect on these systems with great detail. Further- 
more, experiments were the proximity effect on graphene 
could be explored even with atomic scale resolution us- 
ing STM are underway At present only results for 
the total density of states in superconductor-graphene- 
superconductor structures have been presented [lfj ■ 

The present work is aimed to study in detail the inter- 
face between the superconductor and the graphene sheet. 
To this end we shall describe the electronic structure of 
graphene at the level of the tight-binding approximation. 
This description allows us to analyze the superconductor- 



graphene interface more microscopically as compared to 
a description where the continuous limit leading to an 
effective Dirac-Bogoliubov-De Gennes equation is taken 
from the start [ll|. In the continuous description it is 
usually assumed that the presence of the interface do not 
couple different valleys of the graphene band structure, 
which could not be the case in an actual experimental sit- 
uation. Moreover, when the study is focused on finite size 
graphene sheets, a strong dependence on the geometry of 
the edges appears. Thus, different symmetry directions 
will have distinct behavior [12]. For zigzag edges zero- 
energy surface states appear [12j which could hide the 
effects of the coupling to a superconductor. 

In this work we will concentrate on interfaces defined 
along an armchair edge. We propose two different mod- 
els for this interface: the first one assumes that graphene 
is coupled directly to a bulk superconducting electrode 
which does not maintain the honeycomb structure of the 
graphene sheet; the second model studies the possibil- 
ity that one superconducting electrode on top of the 
graphene sheet induces a finite pairing amplitude and 
shifts the Fermi level of the graphene sheet far away from 
the Dirac point. As we discuss below, the two models 
lead to different behavior of the Andreev reflection prob- 
ability as a function of energy, wave vector and doping 
level. We further analyze several aspects of the spectral 
properties of the graphene layer within the two models 
both for the finite and the semi-infinite case. 

The rest of the paper is organized as follows: in Sec. 
II we introduce the tight-binding model for a graphene 
layer and we show the analytic expressions for the Green 
functions for a semi-infinite and a finite layer. In Sec. Ill 
the two different models for the interface with a super- 
conductor are defined and a general expression for the 
self-energy, which provides the basis for the calculations 
of the following sections, is obtained. In Sec. IV we 
study the model dependence of the Andreev reflection 
processes. We also study, in Sec. V, the influence of the 
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FIG. 1: The honeycomb structure of a graphene sheet 
is formed combining two triangular sublattices, denoted A 
(dots) and B (open dots). The unit cell on each horizontal 
line includes one atom from each sublattice. The axis selec- 
tion used in this work is indicated. 



different interface models on the local density of states 
of a finite graphene layer coupled to a superconductor, 
analyzing in particular the minigap which is induced in 
the case of metallic layers. Results for the spatially re- 
solved density of states for a semi-infinite graphene layer 
are presented in Sec. VI. The paper is closed with some 
concluding remarks. 



DESCRIPTION OF ISOLATED GRAPHENE 
LAYERS 

For the description of the electronic states in a de- 
fect free graphene layer we shall adopt the tight-binding 
approximation, i.e. we use a model Hamiltonian of the 
type H = t g E<y>,<7 C £ i CT + e E JCT 4^ where t g de- 
note the hopping element between nearest neighbors car- 
bon atoms on the hexagonal lattice and e is a uniform 
site energy level which allows to vary the level of doping 
(e = corresponds to the undoped case). The dispersion 
relation for the translational invariant case is given by 
E(k,q) = e±^|l + 2e l « a cos(fca/V3)|, where k and q de- 
note the wavevector in the x and y direction respectively 
and a is the lattice parameter defined as shown in Fig. 
[I] (as can be seen a = 3ao/2, where ao is the interatomic 
distance). For the undoped case the Fermi surface col- 
lapse into two nonequivalent points at the Brillouin zone 
corresponding to (k,q) = (±27r/V3a, 0). The dispersion 
relation close to these points can be linearized with a 
slope t g a which fixes the Fermi velocity, vp. 



Green functions for a semi-infinite armchair 
graphene layer 

An essential ingredient for describing the interface be- 
tween a graphene layer and other material is a good de- 
scription of the electronic Green functions at the edges 
of the layer. We concentrate here in the derivation of the 



edge Green function for a semi-infinite graphene layer 
with armchair orientation. We assume that there is 
translational symmetry in the direction parallel to the 
edge (y). The semi-infinite system can be decomposed 
into lines of sites in the y direction which are coupled 
by hopping elements with the neighboring lines on the x 
direction. The unit cell on each line includes two sites 
corresponding to each hexagonal sublattice that are de- 
noted by A and B (see Fig. QJ. These sites are coupled 
by a hopping element t g within the unit cell. Thus, the 
cell Hamiltonian is given by 



h = 



tn 



The hopping elements between neighboring lines cou- 
ple also sites of type A with sites of type B but should 
include a phase factor e qa due to the displacement of 
the cells in the y direction. The hopping matrix in the 
A-B space (both in the forward and in the backward di- 
rection) can be written as t(q) = t g ll(q), where 



U(q) 







The self-similarity of the semi-infinite system with one 
additional line of sites leads to the following implicit 
equation for the edge Green function 



(q,u)= Lui-h-t 2 q U(q)g(q,uj)U(q) 



Hereafter we implicitly assume that u> stands for u>±ir] 
and that the limit r\ — > is taken to obtain the retarded 
or the advanced component respectively. We can now 
define g = Ug which satisfies the simpler equation 



9 = 



X(q,u)—£g 



where X(q, uS) = (u>I — h)U(q). 

To obtain an explicit expression for g it is useful to 
perform a basis rotation in order to diagonalize the ma- 
trix X. The general form of this rotation is R = R1R2 



where Ri 



^iqa/2(7z 



a z being the z-Pauli matrix acting 



on the sublattice space and 



R2 = 



1 



\J— 2 sin a 



-ia/2 



ie- m l 2 ie ia / 2 



where cos a — t g sin qa/ (lj — e) . The eigenvalues of X are 
X\ 2 = — t g cos qa ± \ {oj — e) 2 — t 2 sin 2 qa. We thus get 
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where cos0i,2 = %i,2/(2t g ). The eigenvalues e 4 ^ 1 and 
e -i02 h ave been chosen so that the resulting Green func- 
tions have the proper behavior when the frequency goes 
to infinity. 



Finite graphene layer 

Starting from the results of the previous section one 
can obtain the Green functions of a finite graphene layer 
by introducing a perturbation consisting in breaking the 
bond between the N-th line and its neighbors on the 7V+1 
line. From Dyson's equation we obtain the following set 
of coupled equations 



- F 

9n,u 



9n,n 
9N,n 



9n,N+ng Ni . 
- 9N,N+iig 



jV,n> 



(2) 



where the superindex F stands for the finite system and 
the subindexes i, j indicate the lines within the layer. On 



the other hand the elements <jn,n+i can be expressed as 
9n,n+i = 9N t N^9N+i,N+i- We now use that 51,1 = g 

and 9n,N — 9N,n — (gt) gn,n, where g corresponds to 
the surface Green function for the semi-infinite system 
derived in the previous section. Also we have (Jn+i.n+i = 



[9 - tg^ >N t\ 1 and 9n,n = i - (gt) 

which allows to obtain 



g, 



~F 
9n,n 



*F 

g n ,N 





-1 r 


1 (//')-' 





,2(JV+1) 



1 -1 






i-(gt) 



1 - {gt 

f , -s2(iV+l) 

1 - (gt) 



2(N-n+l) 



(3) 



1 -1 








i- 





(gt) g- 



Making use of the rotation matrix defined in the pre- 
vious section these quantities can be written in the fol- 
lowing rather simple form 



J 



~F 



-F 

g n ,N 




sinn0! sin (N— n+l)0i 
sin 0i sin(A^+l)0i 





1 , sin n<p\ 

/ sin(7V+l)0i 



sin n<p2 sin (N — n+l)4>2 
sin <j>2 sin (7V+1)02 





sin ri(f>2 
sin(Af+l)0 2 



R^U. 



R~ X U 



(4) 



One can have the expression for the borders of the 
layer setting n = 1 or n = N. Then, the eigenvalues of 
the Green functions become sin N fa/ sin (N + 1) fa and 
sin fa/ sin(N + I) fa, with i = 1,2, for the g xl = §n n 
and gf N — gfj x cases, respectively. This expressions are 



equivalent to those for a finite tight-binding chain 15 1. 

The poles of these Green functions determine the spec- 
tral properties of the layer. These poles are fixed by 
the condition sin (N + 1) fa t 2 — 0, which is satisfied 
by 0i,2 = rrnr/(N + 1), where m is an integer. One 
can associate this condition with the quantization of the 
transverse momentum which is used in the continuous 
model for describing armchair nanoribbons (l^ |. At the 
charge neutrality condition the existence of zero energy 



states requires </>i.2 = ±27r/3, which can only be satis- 
fied for N — 3p + 2 (in a more compact notation for 
N mod 3 = 2). Therefore the layers can be classified into 
metallic, for the N mod 3 = 2 case, and insulating for 
the other cases (N mod 3 — 0, 1). In the insulating cases 
the gap in the spectrum is 2E g , where E g ~ ttHvf/3L, 
L = Na/\/3 being the length of the layer. It should be 
noted that electron states in the metallic case are dou- 
bly degenerate, while the degeneracy is removed in the 
insulating cases [13]. 
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MODELING THE 
GRAPHENE-SUPERCONDUCTOR INTERFACE 

One of the aims of the present work is to analyze dif- 
ferent ways to describe the interface between a graphene 
layer and a superconductor. In the recent literature it 
has been assumed that a superconducting electrode de- 
posited on top of graphene induces a finite pairing ampli- 
tude and introduces a finite level of doping which shifts 
the Fermi level of the graphene layer far away from the 
Dirac point 3, 4]. This heavily doped superconduct- 
ing graphene (HDSC) model provides a simple boundary 
condition for the effective Dirac-Bogoliubov-De Gcnncs 
equations describing the interface. 

One can alternatively imagine that the graphene layer 
is coupled directly to a bulk superconducting electrode 
by means of a sharp interface which breaks the coherence 
between the two graphene sublattices. We shall refer to 
this case as the bulk-BCS model. We wish to analyze the 
differences between the two models. 

The presence of superconducting correlations requires 
to introduce the Nambu space describing electron and 
hole propagation within the graphene layer. All Green 
functions acquire a 2 x 2 structure in Nambu space. For 
the uncoupled graphene we have 



9 



9e 

9h 



(5) 



where g e corresponds to the propagators obtained in the 
previous section and gn is obtained from g e by changing 
e — ► — e and t g — > — t g (notice that we use the hat symbol 
to denote the sublattice space while the check symbol 
indicates the Nambu space). 

The effect of the coupling with the superconducting 
electrode can be introduced by means of a self-energy 
S which renormalizes the uncoupled Green functions. 
Thus, the local Green function on the graphene edge at 
the interface is determined by [g^ 1 — 

In the case of the bulk-BCS model the self-energy is 
momentum independent and does not have a structure 
in the sublattice space; i.e. 



£ = Trp s t c f z g B csT z (8 /, 



(6) 



where f z is the z Pauli matrix in Nambu space and t£ is 
the mean square hopping element between the graphene 
layer and the superconducting electrode. This quantity 
controls the value of the parameter — t^Ttp s /t g which 
characterizes the quality of the the interface, p s being 
the density of states. On the other hand, gscs — 9sl + 
f s f x stands for the dimcnsionless BCS Green function, 
i.e. g s = -uj/VA 2 -uj 2 and f s = A/VA 2 -uj 2 . 

We would like now to derive the expression of the self- 
energy within the HDSC model. We first notice that 



in the heavily doped limit one has |e| >> \oj\,t g sin qa 
for the relevant range of frequencies and q values. Thus 
in R 2 one has a — ► tt/2. In the new basis the system 
is equivalent to a tight-binding chain with site energies 
e ± t g and local pairing fixed by A. Although the exact 
Green functions for this system is rather complicated, for 
low energies it can be approximated by 



Tip + 

irp- 
J_f e + t g 
2*2 V e 



gBcs 




(J) 



where irp± = y 1 — ^ ^ /t g . By further taking the 

approximation |e| << t g and transforming back into the 
site representation we obtain the following self-energy 



-T z gBcsT z 



- 1 
- f z 
2 



(8) 



Notice that, in contrast to the first model, the self- 
energy in the HDSC model does exhibit a structure in 
the sublattice space. However, it satisfies the condition 
detE/ig = 1. This structure turns out to be of impor- 
tance in connection to Andreev reflection as discussed in 
the next section. In the following we will use a general 
form of the model self-energies, which can be expressed 

as S = fif z gBCST z ® I — jf z (g) a x . Thus, appropriate 
values for f3 and 7 will correspond to the different mod- 
els (i.e. 7 = with arbitrary (3 for the bulk-BCS model, 
while (3 = a/3/2 and 7 = 1/2 in units oft g , for the HDSC 
model). 



ANDREEV REFLECTION AT A 
GRAPHENE-SUPERCONDUCTOR INTERFACE 



The Andreev reflection is the basic mechanism for the 
conversion of a quasiparticle current into a supercurrent 
at the interface between a normal metal and a supercon- 
ductor. In the case of a graphene-superconductor inter- 
face like the one we have described in the previous section 
there are two channels for the incident electrons with a 
given wave vector q corresponding to the states which 
diagonalize the X matrix. The reflected hole can be in 
either of these two channels. Our microscopic theory can 
thus describe a more general situation than the idealized 
model for the interface used in Q which assumes only 
one channel for the reflected hole for a given wavector. 

The Andreev reflection amplitudes can be expressed in 
terms of Green functions. Generalizing previous works 
1^ . 17 1 we can derive the expression 



f A (q,u}) = 2iA 1 J 2 It 



I-gT, 



A 



1/2 



(9) 



where A e>h (q,u>) = (g e ,h{q, w) - gl ;h (q, w)J /2i Using 
the general form of the model self-energies discussed in 
the previous section allows us to reduce the expression of 
r A to 

f a = 2iiy 2 /3/ s J - /% (g e + g h ) - 7 (^o^ - a x g e ) 



<P 2 +l 2 )9h9e 



A 



1/2 



(10) 



This expression becomes particularly simple when e = 
because g e = g h , f e = ~f h and f' e = -f"\ So f A is 
an scalar quantity given by 



A(3f s (lm 5 2 - |/ - /' 



1-Tr 



(fa 



I T l&x 9e,h 



-(/3 2 + 7 2 ) 



where Tr 



[l3gsi T 
(9 2 -ff)- 



9e,h 



= [W9s9 



V) det ge,/, 

(11) 

j(f + /')] and 



det g e<h 

In general doping conditions (i.e. when e ^ 0) ta is 
not a scalar within the bulk-BCS model. The eigenvalues 
of ta^a Si ye l ne Andreev reflection probability decom- 
posed into two eigenchannels for each wavector q. The 
evolution of these eigenvalues for fixed u> and increasing 
e as a function of q is shown in the left panel of Fig. [5] 
Their maximum value is reached for (3—1 and it never 
exceeds ~ 0.76 at normal incidence. Within the HDSC 
model, however, ta remains scalar for arbitrary doping 
and it always reaches the unitary limit at q = (see right 
panel of Fig.©. 

It is also interesting to analyze the physical charac- 
ter of the Andreev reflection in the two models. The 
information on how the eigenchannels of the uncoupled 
structure are connected by an Andreev process is con- 
tained in the matrix R e rARu- As ta is a scalar within 
the HDSC model, the channel mixing is determined by 
R e R^ 1 . This is also the case for the bulk-BCS model at 
zero doping. In this case we have = 7r — a e and thus 
electrons injected in one channel emerge as holes in the 
opposite one. The momentum in the y direction is con- 
served in this process and therefore this type of reflection 
corresponds to what has been described in 0] as specu- 
lar Andreev reflection. On the other hand for e ^ and 
lo — ► we have ah — a e which corresponds to the usual 
(retro) reflection where holes are reflected on the same 
channel as the incident electron. For intermediate dop- 
ing situations both type of reflection would be present 
although with a dominance of specular (retro) reflection 
for oj > e (u> < e). 

To complete the analysis of the present section we have 
computed the conductance per unit length due to An- 
dreev processes, given by [18( 



Gar = x h I - dqTr ^ {q ' eV) ^ {qi eV °J • " 2 ! 




-0.2 0.2 
v p q/A 



v_q/A 



FIG. 2: Andreev reflection probability on the two eigenchan- 
nels as a function of the parallel momentum q for fixed energy 
lj = 0.2A. The left panel corresponds to the bulk-BCS model 
(with (3 = 1) and the right one to the HDSC model. The 
different curves correspond to different values of the doping 
level e: (full lines), 0.1 (dashed lines), 0.3 (dotted lines) and 
0.4A (dashed-dotted lines). 




eV/A 

FIG. 3: Total conductance per unit length due to Andreev 
reflection Gar normalized to the conductance per unit length 
of a ballistic graphene layer go (V) = ^j-(eV+e)/(irfiVF). The 
results for the HDCS model (full lines) and for the bulk-BCS 
model (dashed lines) are compared with increasing doping 
level e = 0, 0.2, 0.4 and 0.6A. 



The results for Gar within the two models are shown 
in Fig. [3J As suggested in Ref. 0] we normalize the 
result by the conductance per unit length of a ballistic 
graphene sheet, which in the low energy limit is given 
by go(V) — %-(eV + e)/(TrhvF)- As can be observed, 
at the charge neutrality point the HCSD model yields a 
maximum of ratio Gar/ go — 2 at zero voltage, which 
drops to ~ 1.33 for e > 0, in agreement with the results 
of Ref. ||. This ratio is of the order of ~ 1.63 in the 
bulk-BCS model with (3=1 regardless of the doping 
level. The qualitative behavior of the conductance with 
e is similar in both models and agree with the results of 
Ref. (|. 
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PROXIMITY EFFECT ON A FINITE 
GRAPHENE LAYER 

Using the previous results one can analyze the effect 
of the coupling with the superconductor on the spectral 
properties of a graphene layer of finite size. Again, we 
will focus on the differences between the bulk-BCS and 
the HDSC models for the interface. ^From the Dyson's 
equation it is straightforward to obtain the Green func- 
tions at the edge of the layer (labeled as 1) when the 
coupling to the superconductor is introduced on the op- 



J 



posite edge (labeled as N). Then, for an arbitrary line n 
inside the layer we have 



G n ,n — 9n,n + 9 



n,N 



G n ,n — <?n,n + 9n,N I ~ 9n,N^ 



xf 

9N,m 



(13) 



where S stands for the general form of the self-energy 
introduced in Sec. III. We can further reduce this ex- 
pression to 



[/33,-7<^ + (/3 2 +7 2 )&] D e 



[Pgs+i<T x + (P 2 + i 2 )g^ N ] 



9n, 



(14) 



The quantities D e and £>h have the following general 
form, expressed in the graphene subspace 

De = I -09s(9e +9h)~'Y(9eVx-Vx9h) 

-(P 2 + i 2 )g!g F h 

Dh = I - P9s{9e +9h)-l{9h^x-^x9e) 

-(fi 2 +l 2 )9^ 

As in the previous section, in the limit where e = 
these two denominators become equal and simplify to an 



scalar, D, because (3g s (ge + 9h) + l(9e &x 



Tr 



(figslT 7<5"x) 9e,h > and 9e9h = det <? F . Then 



) 



D = 1 - TV [(fig J T l* x ) £ h ] - (/? 2 + 7 2 ) tet£ h 

(15) 

Using the expression of the finite layer Green functions 
given in section II and the rotation given in section I one 
can easily show that 



det[ff e F ,J = - 



sin(7V0i)sin (iV</> 2 



t 2 sin [(AT - 



| sin [(N + l)</> 2 ] 



Tr 



(fig s i T 70"x) 9e,h 



1 



t g sm a 



sin (Nfa) sin [(N + 1)0 2 ] [f3g s - 7 sin {a + q)] 



sin [(N + l)(j>i] sin [(N + l)4> 2 
sin {N(j) 2 ) sin [(N + l)0i] [(3g s + 7 sin (a - q)} 
sin [{N + 1)^] sin [(N + l)<j> 2 ] 

I 



The zeroes of D determine the poles of the coupled sys- 
tem Green functions. From them one can thus analyze 
the distortion of the spectrum due to the superconduct- 



J 



ing proximity effect. For the charge neutrality case, this 
zeroes can be obtained from the expression 



-(3iv 2 



t 2 sin 2 qa^/A 2 



- ur 



+ 



sin [(N + 1)0!] sin [(N + + ((3 2 + 7 2 ) sin (Nfa) sin (Nfo) 
sin (Nfa) sin [(N + l)</> 2 ] - sin (N(j) 2 ) sin [(N + 1)^] 

7 sin (a + q) sin (A^i) sin [(N + l)</> 2 ] + sin (a - q) sin (Nfc) sin [(N + 1)0 X ] 



sin a 



sin (Ncpi) sin [(N + l)<f> 2 ] - sin {N<j) 2 ) sin [(N + 1)^ ] 
I 



(16) 



The spectrum corresponds to a series of subbands 
which disperse quadratically as a function of q in the 
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FIG. 4: Evolution of the lowest energy level Eg of a finite 
graphene layer coupled to a superconductor as a function of 
(3 within the bulk-BCS model for L = 40£ (left panel) and 
with the length of the layer L within both the bulk-BCS (at 
j3 = 1) and the HDSC models (right panel). In the case of 
the bulk-BCS model for arbitrary (3 three different behaviors 
are found depending on the N mod 3 value. On the contrary, 
a universal behavior of E* is found within the HDSC model 
regardless of the N value (see text). 



small q limit. As in the uncoupled case the precise 
form of the dispersion relation depends on the value of 
N mod 3. In the case where the uncoupled layer is metal- 
lic (N mod 3 = 2) the coupling to the superconductor in- 
duces a minigap in the lowest band. The existence of this 
minigap in the spectrum is similar to what is found for 
diffusive conductors and can be associated to the pseudo- 
diffusive behavior of graphene at the charge neutrality 
14| . For the cases N mod 3 = 0, 1 the uncou- 
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point 

pled layer is insulating and the coupling of the supercon- 
ductor just leads to a renormalization of the gap in the 
spectrum. We shall denote by E* the lowest energy level 
for all three cases. 

The dependence E* as a function of N and the inter- 
face parameters (3 and 7 can be obtained from Eq. (fltj|) 
with q = 0. For large N this level decreases as 1/N with 
a prefactor which depends on N mod 3 and the interface 
parameters. Fig. 0] describes the behavior of E* both in 
the bulk-BCS and in the HDSC models. The left panel 
shows the lowest energy state within the bulk-BCS model 
as a function of the interface transparency parameter 
for fixed N. The three cases N mod 3 = 0, 1, 2 are shown. 
As can be observed, in the metallic case the minigap 
evolves from zero at (3 = to a maximum value at j3 = 1 . 
On the other hand, the two insulating cases exhibit dif- 
ferent behavior. While the starting value at (3 = is 
fixed by E g in both cases, in the case N mod 3 = it 
decreases until it reaches the same value as the one of 
the N mod 3 = 2 case for j3 = 1. On the other hand E* 
remains approximately constant for N mod 3 = 1. This 
behavior indicates that the proximity effect is almost neg- 
ligible in this case. The right panel shows the behavior 
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FIG. 5: Spatial variation of the LDOS on a finite layer 
(L — 9£) in the HDSC model. The plots on the top panels 
correspond to the lines inside the layer with n mod 3 = 1,2 
and those on lower ones correspond to n mod 3 = 0. The 
uncoupled layer has a metallic behavior (JV mod 3 = 2). The 
LDOS for this case at the edges of the layer (top left panel) 
and at the center of the layer (lower left panel) is plotted 
with dashed lines. When coupled to the superconductor, the 
LDOS is modified by the appearance of a minigap (denoted 
as mo in the pictures) and with the breaking of the uncoupled 
bands into a pair of subbands (mi and ni2 in the pictures). 
The evolution of the LDOS along the layer is shown in the 
right panels. The results are normalized to the LDOS of a 
bulk graphene layer with zero doping at ui = A, denoted by 
po- 



of the lowest energy state as a function of N both in the 
bulk-BCS model with (5 = 1 and in the HDSC. The re- 
sults are universal (i.e. independent of the ratio A/i s ) 
when plotted as a function of L/£, where £ = hvp/itlS. 
is the superconducting coherence length. It is interest- 
ing to note that while in the bulk-BCS model two lim- 
iting 1/L curves, corresponding to N mod 3 = 0, 2 and 
N mod 3=1, appear, in the HDSC model E* lay on the 
same 1/L curve regardless of N (dashed line in Fig. . 



Local density of states 

We define the electronic local density of states (LDOS) 
on a line n within the graphene layer as 



(27T) 



dqTrlm G n ^ n (q,uj) 



(17) 



which has been normalized to one electron per site and 
spin. The LDOS thus defined is measured in units of 
a/Hvp- However, to study the proximity effect it is more 
convenient to normalize the LDOS with the density of 
a bulk graphene layer with zero doping at w = A, po, 
which for A <C hvp/a is given by A/2ir(a/hvF) 2 ■ The 
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results thus obtained do not depend on the choice of the 
ratio A/tg used in our tight-binding calculations. 

The LDOS on a metallic layer (N mod 3 = 2) is shown 
in Fig. O The results for the coupled case within the 
HDCS model are compared with the results for the un- 
coupled case. It is typically found that the number of 
singularities in the LDOS (associated with the number 
of subbands) in a given energy interval is doubled as 
compared to the uncoupled case. This effect is due to 
the breaking of the double degeneracy of the bands due 
to the coupling with the superconductor. The LDOS 
also exhibits an oscillatory behavior with the position on 
the layer. This behavior reflects the properties of the 
electronic wave functions and, as in the case of isolated 
nanoribbons 12j , is distinct for lines with n mod 3 = 1,2 
and lines with n mod 3 = 0. We thus illustrate these 
cases separately on the top and on the lower panels of 
Fig. El 

The right panels of Fig. [5] show the evolution along the 
layer of the LDOS close to the singularities. We illustrate 
this evolution at three different energies corresponding to 
the lowest first singularities, indicated by m = 0, 1 and 
2 in Fig. [5l For reference we also show the spatial vari- 
ation of the LDOS close to the first singularity for the 
uncoupled case. In this case the LDOS reaches a maxi- 
mum value at the edges of the layer and a minimum in 
the middle for lines n mod 3 = 1,2, while the opposite 
behavior is found for n mod 3 = 0. In the coupled case 
one can still identify the singularities with the oscillation 
pattern in the LDOS but it does no longer reach an ex- 
treme value at the edge of the layer in contact with the 
superconductor. 

As final remark we note that in the case of insulat- 
ing nanoribbons the coupling to the superconductor just 
induces a shift singularities in the spectrum but do not 
change their number. This is due to the nondegenerate 
character of the bands of the uncoupled layer. 



PROXIMITY EFFECT ON A SEMI-INFINITE 
GRAPHENE LAYER 

The results of the previous section can be extended to 
analyze the spectral properties of a semi-infinite graphene 
layer coupled to a superconductor. The local Green func- 
tion on a line n is given by 



fjn 



Cj,, 



- 1 ~t 2 d 

,,n L q-J 



where G n ^ n is the Green function for a finite graphene 
layer coupled to a semi-infinite superconducting layer ob- 
tained in the previous section and g is the Green function 
for the edge of the semi-infinite layer. 

Fig. [S] illustrates the spatial variation of the LDOS on 
a semi-infinite graphene layer and the effect of varying 




FIG. 6: LDOS on lines of type n mod 3 = 1 for a semi-infinite 
graphene layer coupled to a superconductor within the bulk- 
BCS model (lower panels). The upper panels show the cor- 
responding results for the uncoupled case. The plots on the 
left correspond to the undoped case while those on the right 
correspond to e = A. 



the doping level within the bulk-BCS model. For refer- 
ence we show the LDOS for the uncoupled case on the 
upper panels. As can be observed, the uncoupled LDOS 
exhibits long wavelength oscillations on the ~ Hvf/\uj\ 
scale on top of the characteristic V shape behavior. These 
oscillations are a surface effect which decreases in ampli- 
tude inside the layer as shown in Fig. [7J where the LDOS 
profile at cj = 2A is plotted on a larger scale. A similar 
effect has been shown to occur in the case of nanoribbons 
with zigzag edges [l9| . 

The superconducting proximity effect is manifested by 
the appearance of sharp peaks in the LDOS for energies 
|cj| ~ A (lower panels on Fig. [5]). These peaks distort the 
V shape density of states, an effect which decays within 
a few times the coherence length inside the layer. The 
small oscillations on the Hvp/\uj\ scale are reduced as 
compared to the uncoupled case but are still observable 
within the bulk-BCS model (indicated by the full line in 
Fig. [7J). 

The overall behavior of the LDOS within the HDSC 
model is very similar, although in this last case the 
hvF/\u>\ oscillations are further suppressed (dashed line 
in Fig. [JJ. 

The right panels in Fig. [5] illustrate the effect on the 
LDOS of a displacement from the charge neutrality con- 
dition by applying a gate potential (finite e). It is ob- 
served that the V shape is essentially rigidly displaced 
while the peaks induced by the proximity effect remain 
fixed at u> ~ ±A. On the other hand, the wavelength of 
the oscillation pattern is in this case set by Kvf/\uj — e\. 
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L/(27lV F /(0) 

FIG. 7: Oscillation pattern on the LDOS of a semi-infinite 
graphene layer at zero doping within the bulk-BCS model (full 
line), the HDCS model (dashed line) and in the uncoupled 
case (dotted line). The energy is fixed at u = 2A. 



CONCLUSIONS 



slightly smaller in the bulk-BCS model with (3 = 1 than 
in the HDSC model. 

For the semi-infinite case the proximity effect mani- 
fests in the appearance of peaks in the density of states 
for frequencies |w| ~ A. These peaks decay rapidly inside 
the graphene sheet for distances a few times the super- 
conducting coherence length £. On the other hand, the 
LDOS keeps its characteristic V shape for zero doping for 
frequencies < A. We expect that these findings can 
be useful to analyze future STM experiments on graphene 
sheets with superconducting electrodes. 
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We have presented a theoretical analysis of the prox- 
imity effect at a graphene-superconductor interface. For 
this study we have first derived analytical expressions for 
the Green functions on an armchair edge of a semi-infinite 
graphene layer and for a finite layer. 

Two models for describing microscopically the coupling 
to a superconductor have been presented. In the first 
model a bulk superconducting electrode is connected di- 
rectly to the armchair edge of a graphene sheet (bulk- 
BCS model). The honeycomb structure of graphene is 
broken and this is reflected in the different behavior of the 
Andreev reflection probability on the two eigenchannels 
of the graphene sheet as a function of the parallel momen- 
tum q. Only for the case of zero doping both eigenchan- 
nels are equivalent. Within this model one can study the 
effect of varying the normal transparency of the interface 
(through parameter 0) . The Andreev reflection probabil- 
ity at normal incidence never reaches unity within this 
model but has a maximum value of ~ 0.76. In the second 
model it is assumed that the superconducting electrode 
induces a finite order parameter on the graphene regions 
underneath. This model thus maintains the graphene 
sublattice structure and, as we have shown, the Andreev 
reflection amplitude ^(q,^) is a scalar quantity for ar- 
bitrary doping and always reaches the unitary limit for 
normal incidence. 

We have also studied the effect of the coupling to the 
superconductor on the spectral properties of finite and 
semi- infinite graphene layers. For finite layers we have 
obtained a simple expression for the energy spectrum of 
the coupled system which can be easily evaluated numer- 
ically. We have shown that a metallic ribbon develops a 
minigap whose size decreases inversely with the length of 
the layer. This effect can be associated with the pseudo- 
diffusive behavior of graphene. The induced minigap is 
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